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Abstract. - 

We perform molecular-dynamics simulations of a molecular system in supercooled states for 
different values of inertia parameters to provide evidence that the long-time dynamics depends 
only on the equilibrium structure. This observation is consistent with the prediction of the 
mode-coupling theory for the glass transition and with the hypothesis that the potential energy- 
landscape controls the slow dynamics. We also find that dynamical properties at intermediate 
wavenumber depend on the spatial correlation of the molecule's geometrical center. 



Dynamics in normal liquid states deals with orbits in phase space. Binary colhsions and 
vibrations, for which the kinetic energy plays an important role, are elementary ingredients 
building the motion on microscopic time and frequency scales. In supercooled states, the 
system gets trapped in potential energy-landscape pockets for long times. This produces a 
separation of the long-time glassy dynamics — corresponding to the interbasin dynamics, i.e., 
to the motion from one pocket to the other — from the microscopic dynamics. 

In recent years the mode-coupling theory (MCT) for the glass transition has been devel- 
oped pp. One of its essential predictions is the independence of the glass-transition singularity, 
as well as of the long-time glassy dynamics, from the microscopic dynamics. This means that 
the glass-transition singularity and associated slow dynamics are completely determined by 
the statistics of the system's orbits in configuration space rather than in phase space |2 ■ Such 
dynamics, determined solely by the equilibrium structure, is referred to as the structural re- 
laxation. This MCT prediction agrees with the hypothesis that the potential energy-landscape 
controls the slow dynamics [3]. Indeed, an interesting characterization of the MCT critical 
temperature Tc in terms of geometrical properties of the landscape has been proposed 

The issue of the independence of the long-time dynamics from underlying microscopic 
dynamics has been addressed for a binary mixture of Lennard- Jones particles [H] . In Ref. |S] , 
molecular-dynamics (MD) simulations have been performed based on Newtonian as well as 
stochastic dynamics, and the same long-time relaxation was observed. One of the main points 
of this Letter is to explore whether such a property holds for molecular systems. This is 
nontrivial because the dynamics in molecular liquids involves rotational as well as translational 
motions, and there exist more than two time scales characterizing microscopic dynamics. Since 
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these microscopic motions depend on molecule's inertia parameters, we report in this Letter 
MD simulations for systems which differ in these parameters. Let us note that MCT for 
molecular systems - both the molecular theory based on the expansion in generalized spherical 
harmonics (2| and the site-site theory |Hj - predicts the inertia-parameter independence of the 
long-time dynamics, and our work also serves as the first test of such a theoretical prediction. 

Another motivation of the present study concerns the unusual wavenumber dependence of 
characteristic time in molecular systems. A strong correlation in the wavenumber dependence 
is usually observed between the so-called nonergodicity parameters (which specify the strength 
of the so-called a relaxation), the a-relaxation times, and the static structure factor. This 
is found, e.g., in hard-sphere system Lennard- Jones binary mixture jlUj . water and 
silica ^lE], and is in agreement with the MCT prediction. The result found in a MD 
simulation study for a model of orthoterphenyl (OTP) 14 is unusual in that such a correlation 
is violated for wavenumbers about half the position gmax of the first sharp diffraction peak. It 
was suggested in Ref. 14 that the unusual feature is caused by the coupling of the rotational 
motion to the center-of-mass motion. By changing inertia parameters, it shall be examined 
whether it is the coupling to the geometrical center or to the center of mass which matters. 

We consider a model for OTP designed by Lewis and Wahnstrom (LW) and a modified 
model which differs only in the mass distribution. The LW molecule is a rigid isosceles triangle. 
Each of the three sites represents an entire phenyl ring of mass m « 78 amu, and is described 
by a Lennard- Jones sphere. Extensive studies on the dynamics and thermodynamics of the 
original model have recently been performed |14[I16) . In the modified model, the mass of the 
central site, m^, and that of the two adjacent sites, rUs, are chosen so that rric/ms = 16 but 
with the total mass being unchanged, rric + 2ms — 3m. The center-of-mass position in the 
modified molecule nearly coincides with the position of the central site. The thermal velocity 
for the center-of-mass translation is the same for both models. However, the thermal angular 
velocities are increased for the modified model relative to the original one. Thus, one expects 
a considerable speeding up of the reorientational motion, and, via the rotation-translation 
coupling, a considerable change in the complete microscopic dynamics. We have studied, 
performing MD simulations, both the original and modified systems composed of = 343 
molecules at the density p — 1.083 g/cm'^, with the same procedure described in Ref. |16|. 

Three normalized density correlators shall be considered: 

<l>fit) = F^it)/F^{0) for X = N, Z, and Q. (1) 
Here, F^{t) are defined by F^{t) = {p^ (t)* {0)) / N in terms of the linear combinations 

P^~ = ipl+P}+P^/^, p} = {2pl~p}~pl)/V6, p^^ {pi- pi) 1^2, (2) 

of the site-density fluctuations — Yl!i=i exp(i(f • ri) {a = 1, 2, or 3) in which denotes 
the position of site a in ith molecule. (We use the convention that a = 1 refers to the central 
site and a = 2, 3 to the two adjacent sites.) Assuming equal scattering lengths, the correlator 
(/)^(<) is directly related to the cross section as measured in the coherent neutron scattering. 
The functional forms of ^^(t) and 0^(i) have been chosen so that their small- wavenumber 
limits reduce to the Ist-rank reorientational correlators ^7|. Dynamical features of (j)^{t) for 
the original LW OTP have already been discussed in Ref. |14| . 

We show in Fig. QJlhe correlators 0^(t) at T = 260 K for wavenumbers q — 15.2, 19.0 and 
30.5 nm~^. (We notice that the first peak in the static structure factor — F^{0) is located 
at (/max = 14.6 nm~^.) Circles and solid lines denote the results for the original and modified 
systems, respectively, and the time scale of the latter is rescaled by a factor to ~ 0.71, whose 
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Fig. 1 - The correlators (j)^ (t) at T = 260 K for wavenumbers q = 15.2, 19.0, and 30.5 nm . Circles 
and solid lines respectively refer to the results for the original and modified LW OTP, and the time 
scale of the latter is rescaled as t/to with to « 0.71. 

meaning will be discussed below. Except for t < 1 ps, curves for the two systems coincide 
within the statistical errors, i.e., the long-time dynamics is, up to the scale Iq, independent 
of inertia parameters. Figure [3 shows corresponding results for (f>^t) and While the 

short-time dynamics exhibits inertia dependence (see insets of Fig.lzJ, the long-time dynamics 
is identical up to the same scale to. 

The findings in Figs. ^ and |21 verify the following formulation of the glassy dynamics, as 
suggested by MCT 1,2: the microscopic dynamics determines a time scale to so that for 
t^ to any correlator (j)x{t) can be written as 

<t)x{t) ^ Fx{t/to). (3) 

Here, the master function Fx{t) is independent of the microscopic dynamics, and the scale to is 
common to all the correlators. The slow dynamics, as described by Fx{t), is solely determined 
by the potential energy V through the Boltzmann factor e~^/'^^^ (with denoting the 
Boltzmann constant), i.e., by the statistics of orbits in configuration space. The microscopic 
dynamics, which depends on inertia parameters, merely sets the scale to for the exploration of 
the potential energy-landscape. Thus, the factor to used in Figs. ^ and |21 can be understood 
as the ratio to = iJiV^ci of *o's for the original (4) and modified {t]^) systems. The region 
where circles and solid lines coincide is the structural-relaxation part as described by Fx{i). 
The speeding up of the rotational motion for the modified system, demonstrated in the insets 
of Fig. 13 means that the potential energy-landscape can be explored faster. This explains 
why ^0 < 1- 

Figure 13 provides another test of Eq. ©, where we plot 0^ (t) for q ~ 10.2 nm^^ at three 
temperatures close to the MCT critical temperature Tc « 234 K TB . The same to ~ 0.71 has 
been used to rescale the correlators for the modified system. Since the microscopic dynamics 
can also depend on T, a smooth variation of to with T is expected: for example, we checked 
that, by slightly modifying the value of to, an even better agreement in the long-time dynamics 
is obtained for T = 300 K. However, the temperature dependence of to is rather small close 



4 



EUROPHYSICS LETTERS 




t[ps] 

Fig. 2 - The correlators <f>'^{t) for q = 15.2 nm~^ (upper panel) and (f>^(t) for q = 19.0 nm~^ (lower 
panel) at T = 260 K. Circles and solid lines respectively refer to the results for the original and 
modified LW OTP, and the time scale of the the latter is rescaled by to « 0.71. The insets show the 
short-time dynamics, without rescaling the time scale for the modified system, on the linear time axis. 



to Tc as can be seen from Fig. El We note that, within MCT, to should be evaluated at 
T = Tc PEl- 

The dynamics in supercooled states exhibits two-step relaxation as shown in Figs. HEl 
the relaxation toward the plateau, followed by the final relaxation from the plateau to zero 
(the a relaxation). The structural-relaxation region, where circles and solid lines coincide in 
FigS-ES is precisely the one for which MCT is derived and its applicability should be tested. 
We notice from Figs. that the approach toward the plateau, for which MCT predicts an 
asymptotic power-law decay ^ t^°- (0 < a < 0.5) [P, is almost completely masked by the 
microscopic dynamics. This implies that the most faithful tests of MCT should be performed 
just near the plateau and in the a-relaxation regime. A similar conclusion was drawn in 
Ref. UH| in a different context. 

The inset of Fig. 01 shows l/r^ of the a-relaxation times for the original and modified 
systems. In this work, are defined by 0^(t^) = 1/e. MCT predicts an asymptotic formula, 
l/r^ oc {T — Tc)^, close to Tc [l]. The inset shows that the asymptotic formula describes well 
the result with Tc ~ 234 K and 7 w 2.76 as estimated in Ref. ^^l- The exponent 7 is uniquely 
related to the so-called exponent parameter A U. The result indicates that, also for molecular 
systems, Tc and A are independent of inertia parameters, in agreement with the prediction of 
MCT. Notice that this conclusion does not depend on the quality of the fit (i.e., the choice of 
the values for Tc and 7) since, after the rescaling by to, the correlators 0^(i) and hence l/r^ 
for the two systems coincide within the statistical errors. 
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Fig. 3 - The correlators (j>^{t) for q = 10.2 nm at three temperatures T = 230, 260, and 300 K. 
Circles and solid lines respectively refer to the results for the original and modified LW OTP, and the 
time scale of the the latter is rescaled by to « 0.71. The inset exhibits the inverse of the a-relaxation 
times as a function of temperature: l/r^ for the original LW OTP are marked with open circles, 
while those for the modified system (rescaled by to ~ 0.71) with filled triangles. Dashed line denotes 
the fit based on the prediction of MOT, oc (T - T^)^, with Tc « 234K and 7 « 2.76. 



Equation Q is based on the ideal MCT, whose validity is restricted to T > Tc. More 
precisely, Eq. Q is valid also below but according to MCT, Fx(i) for T < Tc would 
describe the arrest at the plateau and does not properly account for the a relaxation in 
the simulation results, e.g., the T = 230 K result shown in Fig. 01 The ideal MCT predicts 
the inertia-parameter independence of the plateau height, and this can be confirmed from the 
figure. However, the T ~ 230 K result also shows that, even for T < Tc, the a relaxation 
seems to be independent of inertia parameters. This implies the possibility of extending the 
formulation JJJ to T < Tc including the whole a-relaxation regime. Such a finding might 
help in developing theories for T < Tc- Let us note that this is also consistent with the 
energy-landscape description whose applicability is not restricted to T > Tc. 

Next we turn our attention to the wavenumber dependence of the a-relaxation time for 
wavenumbers around gmax/2. Figure ^a) exhibits the critical nonergodicity parameters /^'^ 
and the a-relaxation times of the correlators 0^(t) for the original and modified systems. 
Here, f^'^ are the plateau height in the two-step relaxation, and are obtained from the fit based 
on the von-Schweidler law with its first correction as in Ref. Figure^Ja) shows that these 
quantities for the two systems are the same within the statistical errors. One also sees a clear 
correlation in the wavenumber dependence of f^'^ and r^: there exist corresponding peaks 
and minima in these quantities. 

The site-site static structure factor — T'^(O) associated with (j>^{t) is shown in Fig.QJ^b). 
is the same for both the original and modified systems since it depends only on the 
Boltzmann factor. For q > 15 nm^^, the g-dependence of /^"^ and correlates also with that 
of . For g < 15 nm^^, on the other hand, such a correlation is not observed, and there is an 
unusual peak at (7 w 9nm~^ in f^'^ and which does not exist in . This is in contrast to 
results found, e.g., in hard-sphere system [H], Lennard- Jones binary mixture ^OIj water 
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Fig. 4 - (a) The critical nonergodicity parameters f^'^ and the a-relaxation times for the cor- 
relators 4>^{t). The latter are calculated at T = 260 K. Circles and solid lines respectively refer to 
the results for the original and modified LW OTP, and for the latter are rescaled by to « 0.71. 
(b) The static structure factors at T = 260 K. Solid and dashed lines respectively denote the site- 
site static structure factor and the geometrical-center static structure factor S^'^, which are the 
same for both the original and modified LW OTP. For the original system, S^'^ coincides with the 
center-of-mass static structure factor S^'^^. Sq'~'^ for the modified system is denoted as dotted line. 

and silica |12lll3j . for which such a correlation holds for the whole wavenumber regime. 

In Ref. it was suggested that the unusual feature around q~ 9 nm^^ can be interpreted 
as being due to the coupling to the center-of-mass (COM) dynamics. This follows from the 
fact that the COM static structure factor, S^°^ = (l/N) ^. (e-^^ ^*^""""^^"")) with rv^OM 
denoting the COM position, has a peak around q w 9nm~^ as shown in Fig. Elb). For 
the original model, the COM position is identical to the geometrical center (GC) defined 
by f;GC ^ {l/3)J2l=iK''- The latter is a structural property, i.e., independent of inertia 
parameters, and the corresponding structure factor S^'~^ , defined in terms of the GC positions, 
is the same for the two systems in contrast to S^'~'^ {cf. Fig. ^b)). Since Fig. 2fa) shows 
that the peaks in /^"^ and are also a structural property, we find — from the correlation 
between f^'^, , and 5*^^ — that it is more appropriate to interpret the unusual peak around 
q « 9nm~^ as being caused by the coupling to the GC dynamics. Thus, the unusual peak 
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reflects a property of the potential energy-landscape, a subtle interplay of the translational and 
rotational dimensions in configuration space. We notice that such a peak in the nonergodicity 
parameters around q « 9nm~^ can be found in the coherent neutron-scattering result for real 
OTP molecules, where no clear peak shows up in the static structure factor (19^ . 

The position q w 9nm~^ where the unusual peak we discussed occurs is compatible with 
the inverse of the van der Waals radius rw = 0.37 nm for OTP molecule '50", i.e., it is connected 
to the overall size of the molecule. It is interesting to note that a similar unusual peak was 
found in a model for polymer around wavenumbers close to the inverse radius of gyration |21| . 
In Ref. j21| . the unusual peak is interpreted in terms of polymer specific properties. Such 
an interpretation might not be suitable, if the unusual peak - found in OTP and polymer - 
reflects the same physics. Thus, it will be a theoretical challenge to reproduce and explain 
such an unusual feature observed in molecular and polymer systems. 

* * * 
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